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Using detailed balance and scaling properties of integrals that appear in the Coulomb gas refor- 
mulation of quantum impurity problems, we establish exact relations between the nonequilibrium 
transfer rates of the boundary sine-Gordon and the anisotropic Kondo model at zero tempera- 
ture. Combining these results with findings from the thermodynamic Bethe ansatz, we derive exact 
closed form expressions for the transfer rate in the biased spin-boson model in the scaling limit. 
They illustrate how the crossover from weak to strong tunneling takes place. Using a conjectured 
correspondence between the transfer and the decoherence rate, we also determine the exact lower 
bound for damping of the coherent oscillation as a function of bias and dissipation strength in this 
paradigmic model for NMR and superposition of macroscopically distinct states (qubits). 
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I. INTRODUCTION 

Quantum impurity problems (QIPs) have attracted a 
great deal of interest recently. This is because the under- 
lying physics is non-trivial and the models are manage- 
able technically despite their essentially nonperturbative 
nature. In addition, they have a multitude of experi- 
mental applications, including the Kondo effect, quan- 
tum dots, dissipative quantum mechanics, tunneling in 
quantum wires and fractional quantum Hall devices i 

There have been discovered various relations between 
thermodynamic quantities of the anisotropic Kondo 
model (AKM) and the boundary sine-Gordon (BSG) 
field theory models Each of these integrable models is 
of considerable interest and it is remarkable that they 
are closely related. They are both boundary integrable 
quantum field theories with a quantum-group spin at the 
boundary which takes values in standard or cyclic repre- 
sentations of the quantum group SU(2) q & 

The AKM is related to the spin-boson (SB) models 
The BSG model is equivalent to the Schmid model, which 
describes a damped particle in a tilted cosine potential 5 
and has a great variety of applications to charge and 
particle transport including Josephson junction dynam- 
icsi The field theory limit in the AKM and BSG model 
corresponds to the Ohmic scaling limit in the SB and 
Schmid model. Here we are interested in the strong- 
backscattering or tight-binding (TB) representations of 
these models. The equivalence and difference of the vari- 
ous models can most easily be seen in the Anderson- Yuval 
Coulomb gas representation for the partition function, 
which is the grandcanonical sum of a one-dimensional 
classical gas of positive and negative unit charges with 
" log-sine" interactions and overall neutrality. The AKM 
or SB model differs from the BSG or Schmid model by 
a different charge ordering. Principally in the sequel, we 
concentrate on the Schmid and SB model and use the 
language of dissipative quantum mechanics. With use of 
the correspondence relations, our results directly apply 
also to the BSG and the anistropic Kondo model. 



Here we study nonequilibrium quantum dynamics of 
the models mentioned with the rigorous Keldysh ap- 
proach. Upon uncovering functional relations between 
perturbative Coulomb integrals, we discover that at zero 
temperature for every particular transfer rate only paths 
with minimal number of tunneling moves contribute. 
This allows us then to find the characteristic function 
encapsulating all statistical fluctuations of the transport 
process. The results for the Schmid model corroborate 
the findings from the thermodynamic Bethe ansatz. We 
discover, as a second result of current interest, an in- 
triguing functional relation between transfer rates in the 
Schmid and partial transfer rates in the SB model. The 
major result is given in Eq. (J73J). Using this, we derive a 
concise contour integral expression for the transfer rate in 
the SB model [cf. Eq. (|86[) ] , from where we find both the 
weak- and the strong-tunneling expansion for this rate. 
We also make and check a conjecture which relates the 
transfer rate to the decoherence rate in the SB model. 
The ensuing asymptotic series, Eq. H1UU|) . provides an 
exact lower bound for decoherence in the scaling limit as 
a function of bias and damping strength. 

In the following section, we briefly review a general 
Poissonian transport model with transfer rates describ- 
ing direct forward and backward transitions by m states, 
m = 1, 2, 3, • • • . We solve the master equation and 
determine the evolution of the probability distribution. 
In Sec. IIIII we study the Schmid model in the real- 
time Coulomb gas representation, which is based on the 
nonequilibrium Keldysh method. Upon performing a 
cluster decomposition, we show that the BSG dynamics 
at long times represents a realization of the Poissonian 
dynamics studied in Sec.|Hj We derive exact formal ex- 
pressions for the quantum transfer rates in which friction 
and noise effects are clearly separated. The correspond- 
ing investigation for the AKM or SB model is presented 
in Sec. IIVI The weak-tunneling expansion for the trans- 
fer rates at zero temperature are studied both for the 
Schmid and SB model in Sec. Using detailed balance 
relations, and scaling properties holding for a subgroup of 
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the Coulomb integrals, we uncover a multitude of linear 
functional relations between the occuring Coulomb inte- 
grals of same order. Using these relations, we discover 
formidable cancellations in the path sum for the individ- 
ual rates of the Schmid model such that only the paths 
with minimal number of tunneling moves contribute. 
This peculiarity allows us to determine all statistical fluc- 
tuations once the current is known. The corresponding 
discussion is given in Subsec. IV El We also discover an 
intriguing functional relation between the rates in the 
Schmid model and the partial rates in the SB model. 
The findings in Sec. El can be used to shed light on the 
entire regime from weak to strong tunneling, both in the 
Schmid and SB model. In Sec. ED we derive from the per- 
turbative series an exact integral representation for the 
transfer rate in the SB model at zero temperature valid 
for general coupling, bias and damping strengths. From 
that, the asymptotic (strong-tunneling) expansion for the 
transfer rate is derived. In Subsec. IVI CI we propose a 
conjecture about a simple close relation between the re- 
laxation and decoherence rate in the SB model. The con- 
jectured formula for the decoherence rate is in agreement 
with all known special cases and gives a lower bound on 
decoherence for general bias and damping strength. Fi- 
nally, in Subsec. lVIDI we give analytic expressions for the 
leading thermal corrections. 



II. POISSONIAN TRANSPORT AND NOISE 

Consider a general model describing transport of mass 
or charge between sites n = 0, ±1, ±2, • • • via direct for- 
ward and backward transitions by £ sites, £ = 1, 2, • • • . 
The respective weights per unit time (transfer "rates" ) 
are denoted by 7^. Assuming statistically independent 
transitions, the dynamics of the population probability 
P n (t) of site n is governed by the master equation 
00 

Pn{t) = ^ biPn-i(t)+7iPn+e(t) ~ (7/ +7i")-P«(*)] ■ 
1=1 

The Fourier transform P(k, t) = J2 n e * fc ™ Pnif) with ini- 
tial state P n (0) = S n is found from this equation as 



P(k,t)= f[exp [ t (e lkn -l)j+ + t 



e~ ikn - 1) % 



The moments of the probability distribution follow from 
the characteristic function P(k, t) by differentiation, 



*(*)>= 5>"*P n (t) = (-%—) P(k,t) 



d 



(i) 



k=0 



The resulting expressions can be written in terms of the 
irreducible moments or cumulants (n m (t)) c as 

(»(«)) = («(*)>c . 
(n 2 (t)) = (n(t)) 2 c + (n 2 (t)) c , (2) 

(n 3 (t)) = (n(t)f c +3(n 2 (t)) c (n(t)) c + (n 3 (t)) c . 



The cumulants expansion 

00 

lnP(M)= Yl 

m— 1 

leads us to the expressions 



(ik) T 



(n m (t)) c 



(3) 



(n m (t)) c 



fc=0 



*X> m bn+(-l) ro 7n] 



(4) 



71=1 



The function P(k,t) may be conveniently rewritten in 
terms of partial forward/backward currents I„ — 17^ as 



P(k,t) = l[Z+(k,t)Z-(k,t), 



71=1 
OO 



knt (tin / n Y 
£\ 



-tl±/n 



(5) 



exp[i(e 4fc " - l)/±/n 



The characteristic function P(k, t) encapsulates all sta- 
tistical fluctuations of the transport process. To elucidate 
the physical meaning of the above expressions, suppose 
now that the mass transferred per unit time in forward 
direction were the result of a Poisson process for parti- 
cles of unit mass propagating via uncorrelated nearest- 
neighbour forward steps, contributing a current = 7^, 
plus a Poisson process of uncorrelated forward moves via 
next-to-nearest-neighbour transitions contributing a cur- 
rent I2 = 1 etc. Further suppose that the total back- 
ward current were the result of independent Poisson pro- 
cesses with corresponding partial backward currents I~ , 
n = 1, 2, • • • . The final form of the characteristic function 
would then be the expression JSJ. 

With regard to applications of this model to QIPs (see 
below), there is a different interpretation. Expression 
|J5J may also represent independent Poisson processes for 
particles of charge one going across the impurity's bar- 
rier in forward/backward direction, contributing a cur- 
rent If' , plus a Poisson process for particles of charge 
two (or for a pair of particles of charge one) contributing 
a forward/backward current J^, etc. 

Next, suppose that the particle is coupled to a thermal 
reservoir at temperature T . Then the forward and back- 
ward transfer weights are related by detailed balance. 
Assuming that the potential drop per site interval in for- 
ward direction is e > (we use units where H = k B = 1), 
we then have 

7" = e - e/T 7+ ■ (6) 

For classical Poisson processes, all of the transfer 
weights 7„, n = 1, • • • , must be positive. When the 
Poisson processes come about quantum-mechanically, the 
transfer weights may be partly negative (see below). 
This does not spoil conservation of probability, since 
J2„ P(t) = by construction of the master equation re- 
gardless of the forms chosen for the set {7^}- 
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III. NONEQUILIBRIUM QUANTUM 
TRANSPORT IN THE BSG OR SCHMID MODEL 

The noncquilibrium transport through a local 
backscattering potential embedded in a Luttinger liquid 
environment gives a fingerprint of the non-Fermi-liquid 
state. This situation is realized e.g. in resonant tun- 
neling experiments through a point contact in quantum 
Hall devices*^! In the field theory limit, the tunneling 
problem is described by the BSG model, which is a har- 
monic bosonic field on a half-line, and the only interac- 
tions take place on the boundary^ The Hamiltonian for 
strong-backscattering is H BSG = H + H Bound , where 



-in 



^d*[H 2 + (^0) 2 ]+— ^ dxd x 



H 



Bound 



A s cos[<K0); 



(7) 



Here, H represents the bulk and an applied voltage V, 



and H, 



Bo 



is the interaction on the boundary. 



The Schmid model describes a quantum Brownian par- 
ticle moving in a tilted cosine potential^ In the TB repre- 
sentation, the Hamiltonian is (except for a counterterm) 



H s = — ^E + h - c -) ~ e E na ™ a « 

n n 
- E C <* X a E na «°« + H bzth( { x a\ ) • 



(8) 



The first line gives the TB system. Here, A s is the cou- 
pling energy and e is the potential drop or bias energy 
between neighbouring sites (lattice constant a = 1). The 
second line describes the system-bath coupling and the 
harmonic bath. In the Ohmic scaling limit, the spectral 
density J(u>) = 7rE a <K w - Lo a )c 2 a / {2m a L0 a ) is strictly 
linear in ui, J{u) — 2Ku>. The Kondo parameter K is a 
dimensionless Ohmic damping strength, 8 and is the in- 
verse of the Luttinger parameter g or the filling fraction 
v in fractional quantum Hall systems. 

The effects of the Luttinger liquid bulk are in the two- 
point function on the boundary, (e s< ^ 0,t ) e -l ^°'°') T = 
e -0(t). The function Q(t) = ( [0(0,0) - 0(0,*)] 0(0, 0)) T 
directly corresponds to the thermal correlator in the 
Schmid model. Introducing a high-frequency cutoff u> c , 
we get in the field theory or scaling limit 

Q{t) = 2!Hn[(w c /7rT)sinh(7rT|i|)] + inKsga(t) . (9) 

The BSG model is in correspondence with the Ohmic 
Schmid model if we put K = 1/g and e = 2ttV. The 
system (JJJ or (JSJ represents a quantum-mechanical real- 
ization of the general model discussed in Section II. 

The nonequilibrium transport may be found upon em- 
ploying the Feynman- Vernon or Keldysh formalism to 
the calculation of the reduced density matrix (RDM). 
It is convenient to parametrize the sudden moves of the 
system along the forward path q and backward path q' 



in terms of charges {uj — ±1} and {vi = ±1}, respec- 
tively. It is then straightforward to derive the exact se- 
ries expansion in the number of tunneling transitions for 
the characteristic function P(k,t). This series has formal 
similarity to the partition function of charges in a one- 
dimensional fixed volume. Taking the Laplace transform 
P(fc, A) formally corresponds to switching to the analog 
of an isobaric ensemble of charges. Upon tracing out the 
thermal bath, the charges are interacting with each other 
through the complex time- nonlocal interaction Q(t). 

Consider an individual double path or charge sequence 
with n + 2£ time-ordered u charges on the q path at lo- 
cations f i < *2 • • • < t n+ 2i and n + 2m time-ordered v 
charges on the q' path at locations t[ <t' 2 ---< t' n , 2m , 



,±n _ 



{ui,u 2 • • • , u n+2 e; Vi,V2, ■ ■ ■ , v n+ 

2rn / ±n , 



where the total charges satisfy the constraints 

n+2t n+2m 

J2 u j = E v ^ ±n - 

3=1 1=1 



(11) 



The tunneling amplitude factor for this arrangement is 

(12) 



.>■.,::.! = I*) (=± 



and the bias phase factor reads B[a 



±n i _ itp[a e ™; 



with 



(13) 



The interactions for the arrangement (|1L)|) are described 
by the Gaussian influence function 

F[<& =^[lt n }^iW m n ]^[qt n ,qt n } , (14) 

where T\ and describe the selfinteractions of the q and 
q' path, respectively, and J- 2 the mutual interactions, 

, n+2t >. 

Fi[qf n } = expj Y UjUiQ(tj -*»)> , 

^ 3>i=l ' 

{n+2m >. 
^2 VjViQ*^ -*•)}, (15) 
3>i=l ' 

{n+2m n+2l n 
-EE '■■"■Wj '■>} ■ 
3=1 i=l > 

The Poissonian dynamics (JSJ) appears at times * large 
compared to the transition times between the sites. It is 
found upon performing a cluster decomposition of the in- 
dividual charge arrangements contributing to P n (*) . The 
clusters represent the elementary processes of the tunnel- 
ing dynamics. They are charge sequences corresponding 
to irreducible path sections which start and end in diag- 
onal states of the RDM. By definition, clusters are non- 
interacting. Therefore in Laplace space, a succession of 
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clusters factorizes and every interval in a diagonal state 
separating neighbouring clusters generates an additional 
factor of 1/A. Path segments which have interim visits of 
diagonal states can be divided up into a reducible part, 
which breaks down into a succession of clusters of lower 
order, and an irreducible remnant (see below). 

When the Laplace variable A goes to zero, the clus- 
ters become independent of A. As a result, the dynamics 
turns into Markovian behaviour at long times. In this 
limit, the totality of irreducible clusters describing the 
elementary processes — > ±n and 0^0 represent the 
weight per unit time 7^ for direct transitions from arbi- 
trary site m to site m ± n, and the weight per unit time 
7 to stay at the same site, respectively. Since the sys- 



tem with initial state P n (0) 



J n,0 



must finally take real 



occupation of site n, there is an excess of forward moves 
over the backward moves for n > 0. Thus we find from 
the Coulomb gas or path sum (we put S(m, n) = 5 m n ) 



OO I OO OO / -L. 



P M = E E 



i=l ( j i= Q e.=o Jl ' l ' 



(16) 



The transfer "rate" 7+ can be divided up into 7 + (+) 
and 7„ (— ), in which the last charge is positive and neg- 
ative, respectively, 7+ = 7 + (+) + 7«( — )• Next, assume 
that we replace the last charge in 7^(— ), which is either 
ti f = —1 or v { = —1, by the charge v { = +1 or u { + 1. 
Under this substitution, the bias phase and the influ- 
ence factor are left unchanged while the amplitude |Q 
changes sign, and the constraint n in Eq. (fill) changes 
into ri + l, as we see from Eqs. (|11|) - <|15[) . Therefore 
we have 7+H = -7,t+i(+)- Similarly, we get 7~(+) = 

-7n+l( - )- ThuS WC nnd 7n = 7n (±) " 7n+l(±) and 

7 = — 7i~(+) — ), and therefore finally 



00 

7o + E(7,l +7,7) =0 



(17) 



This relation ensures conservation of probability of the 
expression l|16[l. ^2 n P n {t) = 1. Furthermore, the char- 
acteristic function P(k,t) is found to take the form (J5J). 

The irreducible contribution of the charge arrangement 
(jTU|l with excess charge ±n in Eq. ljTT|) gives a partial 
contribution to the transfer "rate" 7„ . We find 



(-1) 



£+m 



A. 



2(n+l+m) 



At 



(18) 



~ n+2l-l n+2m-l 

j n ^ n M^^^Hafa], 



j=i z=i 



where is the irreducible part of the influence func- 
tion. The integrations are over the charge intervals 



Pj = tj + i — tj and p\ = t' i+1 — t\. The additional un- 
restricted integration over the interval t = t' x — t n +2i 
introduces iV„ Am = (j+^^n+Tm)! different possibilities 
of mixing up the time-ordered moves in q^ n with those in 
q'^j™. In general, the partial transfer weight (|18fl is com- 
plex. The complex conjugate counterpart comes from the 
set <|1(J[1 in which the u and v charges are interchanged. 

For given n, the term with £ = rn = in Eq. i|18fl has 
minimal number of tunneling moves. The contributions 
with t, m > have additional pairs of forward/backward 
moves both along the path q and the path q' . The grand- 
canonical ensemble of charge arrangements satisfying the 
constraint (|11|) yields the transfer rate 7^, 



it 



1=0 m=0 „±« 



(19) 



The sum ±n is over the M n -t m — ^ n+ l e ?l ("+ 2 ") 1 

*-~'OL t n,t,iu [n+£)l£l (n+m)\m\ 

different possibilities to order the u charges for fixed n 
and £ and the v charges for fixed n and m. 

The partial transfer weight (|18[) takes up the tunnel- 
ing contributions to the incoherent rate 7^, in which (i) 
the succession of moves in the q and q' path is kept fixed 
separately, and (ii) all possible arrangements relative to 
each other are contained, the latter being due to the unre- 
stricted r-integration in Eq. (|18J) . It describes a coherent 
tunneling process from site to site ±n and takes into 
account M„ ; ^ m different possibilities of visiting virtually 
2(n + £ + m) — 1 intermediate states of the RDM. 

It is straightforward to see that the partial "rate" 
(I18JI obeys detailed balance already. To this, observe 
that the expression J^J is analytic for complex time z 
in the strip > Imz > — 1/T and has the property 
Q(t - i/T) = Q*(t) = Q{-t). Thus we may shift the 
flip times of the backward path as — > t£ — if, where 
< f < 1/T. This results in a shift of the r integration 
path in l|18fl . T — > r — if . Since this path is unrestricted, 
the edge contours perpendicular to the real-time axis do 
not contribute. Hence the integral does not depend on 
the shift f. Putting f = 1/T, the bias factor l|13|) in 
Eq. i|18|) then generates a factor e ±ne / T and the interac- 
tion term T2 transforms into T^i while T-y and T* are 
left unchanged. Next observe that under combination of 
time reversal and charge conjugation both the bias factor 
e lip and the interaction factors T\ and T\ are invariant, 
while T2 transforms back into J^2- Thus we find 



7n 



(20) 



where 57™. is the charge sequence resulting from (f 1 01) by 
charge conjugation and time reversal. 

In the above we have shown that the {u, v} representa- 
tion is useful to identify partial sums of tunneling paths 
for which detailed balance holds already. On the other 
hand, the disadvantage of this representation is that it 
mixes up friction and noise. These effects can be sepa- 
rated clearly in the Ohmic scaling limit using a different 
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representation. To this end, we now introduce charges 
rjj = ±1 describing forward/backward moves along the 
quasiclassical path q + q' and charges & = ±1 represent- 
ing sudden moves along the quantum fluctuation path 
q — q'. We then have the correspondence 



it 

v 



±1 
±1 



{77, £} = {±1, ±1}, 
{77, f} = {±l,Tl}. 



The cumulative charges 



fe=l 



and 



E4 (21) 



fc=i 



measure propagation of the system after j moves in for- 
ward/backward and off-diagonal direction, respectively. 

An individual tunneling path with 2£ moves contribut- 
ing to the rate 7^, where I > n is then parametrized by 
a set of {77, £} pairs at locations t\ < t% < • • • < t^i, 



±n 



{771, Ci; »72,6; • • • ; }±« 



where the charge arrangements are constrained as 

2t 21 



(22) 



(23) 



Next, we introduce dimensionless temperature t? = T/e 
and dimensionless charge intervals T i — t — U) where 
i = 1, • • • , 21 — 1. Then the irreducible influence factor 
for the charge sequence l|22[) may be written asm 



(24) 



The influence phase <^ (77, £ ) adds up the imaginary parts 
of the bath correlations and describes friction, 



21-1 



(25) 



Observe that <f>i depends only on the order of the charges, 
but not on the individual charge intervals r l7 • • • , T 2e _ 1 . 

The amplitude Gg (r, £ ) includes the quantum fluc- 
tuations. The term Gf uhtl \f, £) subtracted from 



G[ c) {T 1 i) = G t (T,i)-G 



(subtr) / -. 



(r,0 



(26) 



is to eliminate the reducible parts of G t (r, £ ), if any. 
Introducing the distance of charge pair {j, «}, 



fe > 



the amplitude takes the form 



21 



Gt(r,£)=exp ( ln 
j>i=i 



sinh(7n?r- i ) 



(27) 



Subtractions are necessary when a cluster splits up into 
neutral subclusters, i.e. when one or several of the cu- 
mulative p-charges are zero. By definition, in the sub- 
tracted term G| Subt1 '' (r, £) all interactions across an in- 
terval associated with zero cumulative charge are disre- 
garded. Hence, Gf ubtI ^ (r, f ) factorizes into a product of 
amplitudes of lower order representing the subclusters. 
Putting, e.g., 1=2 and £2 = —£,1, i-e. p 2 = 0, we have 

2K / a \ 2if 



G 



(subtr 



(r,0 



7n9 



sinh(7T7?r 1 ^ 



sinh(7r7?r 3 ) 



The irreducible amplitude C?» (r, £ ) encapsulates the in- 
teractions of the £ charges and acts as a Gaussian filter 
controlling the quantum fluctuations. 

The residual factor {e/u c ) 2Ke in Eq. QZQ is an adia- 
batic Franck-Condon factor leading to renormalization of 
the tunneling amplitude, ft is convenient to absorb A, e 
and oj c into a single dimensionless parameter, 

K A 

■'s = I — I — ■ (28) 



The full rate 7^ can now be written as 



(29) 



U tt = X [cos^C,±sin^5 £ ] . (30) 

The prime {•••}' is to remind us that the double sum 
underlies the constraints 1)23(1 . There are ^f^- ^ + ^y l ^_ n y 

different charge sequences contributing to Uz: «■ All quan- 
tum fluctuations are in the 21 — f-fold irreducible noise 
integrals (f ) = C £ (£ ) ± iS^f), 

/■OO 

J±(0 = / d 2 ^V G< c) (f, i) e ±W.O , (31) 



where J °° d 2£ 1 r . . . = J °° dr 2£1 • • ■ di^ .... In terms 
of the cumulative p charges, the bias phase l|13|l reads 



2£-l 



e(r,f)= £>. 



(32) 



The integrand of the noise integral C £ (£>^) depends on the 
cosine (sine) of the bias phase and the interactions of the 
"fluctuation" charges are in the factor G* c) (r,£). With 
the subtractions in Eq. I|26|) . all integrals in Eq. (|31|l are 
convergent at large intervals between the charges. 

The functions U^r. are linear combinations of the noise 

integrals C t {^) and S e (£). As a consequence, the multi- 
tude of detailed balance relations (|2C)|> for partial rates 
directly results in a wealth of relations between noise in- 
tegrals jf~ (£ ) of the same order. The coefficients of these 
relations include all the friction effects. We conclude this 
section with the remark that it is only the Ohmic scaling 
limit in which complete separation of friction and quan- 
tum or classical noise takes place. 
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IV. INCOHERENT TUNNELING IN THE 
SPIN-BOSON MODEL 

The other quantum impurity model of interest is the 
one-channel anisotropic Kondo models of spin ^ , which 
describes scattering of fermions at an impurity with in- 
dependent exchange constants pJ« and pJ ± . Here pJ^ 
conserves the polarization of the spin, while pJ ± causes 
spin-flip scattering. The AKM can also be expressed in 
bosonized form, H = H + H spin , where H is given in 
Eq. " 



1 ■"spin' 

JjJ. The boundary Hamiltonian is 



l_l 

spin 



A 



sl; 



e 

2 2 



(33) 



The additional bias term describes a local magnetic field 
energy. In general, the relations between the parameters 
g and A SB of the boundary spin model with the param- 
eters pJn and pJ ± of the AKM depend on the particular 
rcgularization prescription chosen for the singular scat- 
tering potential. iiS 

On the other hand, the boundary spin model l|33() with 
(|7| is the polaron-transformed equivalent of the spin- 
boson (SB) model in the Ohmic scaling limitfi^ 



#SB - 



i SB 



e 1 

2 z 2 



2 / J 



bath 



(34) 



In view of application of the SB model to a damped par- 
ticle in a double well in the two-state limit, we choose 
the basis to be formed by the localized states |R) (right) 
and |L) (left), which are eigenstates of the Pauli matrix 
a z with eigenvalues +1 and —1. The potential drop from 
\L) to \R) is e, and A SB is the tunneling coupling. Again, 
we have the correspondence K — l/g and e = 27rl^. 

A valid dynamical quantity of interest is the expecta- 
tion value of a z at time t describing the difference of the 
populations of the two states, (cr z {t)} = PH.(t) — Ph(t)- 
The dynamics of (a z (t)} may be described in terms of the 
exact general master equation (GME) 



dt 



dt' 



K ia) {t-t') + K (s) (t-t')(a z (t')) 



or equivalcntly in Laplace space 

A(<7 2 (A)) = -k^(X)-k^(X)(a z (X)) , 



(35) 



where the kernels K^ s /^{t) and K^ s / a \X) are even/odd 
under bias inversion. 

In this section we are interested in the regime where the 
kernels decay fast compared to the relevant time scales 
of the relaxation dynamics of the particle. On this as- 
sumption, the GME becomes local in time, 



where 



dt 



7± 



7" - 7 - (7 + +r)(^W) 



kW(t) ± if (a) (r) 



k {s \x = o)±k^ a \x = o) 



(36) 



(37) 



is the forward/backward incoherent tunnneling rate. The 
solution of Eq. H36[l is 



Mt)) = (cr,) eq +[( O -,(0))-<cr,) Bq ] 
where (<7 z ) eq is the equilibrium value, 



-(r+7 )t 



'2/eq 



— — — = tanh ( -—i 



7 T + 7 



\2T) 



(38) 



(39) 



The second form follows from detailed balance. 

The expression <|38|) follows again from a cluster ex- 
pansion for {<7 Z (X)). The obvious difference between the 
Schmid and the SB model is that in the former the time- 
ordered u and v charges may have arbitrary charge order, 
while in the latter they must alternate in sign. In the SB 
model, the 2£+l u charges and 2m + 1 v charges in 1|10|) 
are therefore separately ordered as 



^ = {±, T ,±,--- ,±; ±,t,±,--- ,±} 



(40) 



Thus, the series for the total tunneling rate in {u, v} 
representation is 



oo oo 

e=o m=0 
„±r yl±l 



(41) 



The rate contribution 7^ [.A/vL] is given in l|18fl . where 
A gB takes the place of A s . Thus, the full rate in the SB 
model corresponds to a particular partial rate contribut- 
ing to the rate 7^ in the Schmid model. 

Next, we switch to the {??,£} representation. For the 
set of paths specified by (|4*0"|) the sequence of 77's and £'s 
are restricted such that the cumulative charges r ■ and p- 
(j = 1, • • • , 2 + 21 + 2m) obey the constraints 



< r j < 2 , 
> r i > -2 
Pi = 0, 



for 
for 
for 



a: 



A. m . (42) 
j even . 



In analogy with Eq. (|28|l we introduce the dimensionless 
dressed tunneling coupling as 



''SB 



< \ K A S B 



Formally we may then write 



00 



(43) 



(44) 



where Wf' — , and where the double sum in 13U|) 

underlies the SB constraints 1)42(1 . The r\ sum in U^ t is 
easily performed, yielding 



Wf = 



cosV^) ^ ± 6 tan^-ftf)^] , (45) 
fc}" 



where {•}" reminds us of the SB constraint £ 2fc — — £ 2 fc-i 
(k = !,•■•, t) [cf. Eq. 03 . 
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V. WEAK-TUNNELING EXPANSION AT T = 



B. The order A 



In this section we investigate the weak-tunneling ex- 
pansion for the tunneling rates in the Schmid and SB 
model. We restrict the attention to the scaling limit and 
zero temperature. Then the charge interaction is loga- 
rithmic at all distances, Q(t) = Q {t) where 



) (t) = 2K\tl(t) + i ttK sgn(r) . 



(46) 



Consequently, the noise integrals (£) given in Eq. ijol"|) 
diverge at short distances between neighbouring attrac- 
tive charge pairs for K > ^ . Now we first focus the atten- 
tion on the case K < ~ , in which all noise integrals are fi- 
nite. The way to regularize the integrals for K > i with- 
out introducing nonuniversal cutoff-dependent physics 
shall be discussed at the end of Subsection IVI Al 

At T = 0, backward transitions to higher wells are 
absent, i.e., 



In [Ol-l, 







for all 



(47) 



as follows from Eq. (|2U|) . In the sequel, the u charges 
are depicted as + and — , and the v charges as © and 
G. Complex conjugation is synomynous with exchange 
of the u and v charges, and charge conjugation means 
change of the signs of the charges. 



In order Ag in the Schmid model, there is (i) the pos- 
sibility for direct tunneling transitions to the next-to- 
nearest-neighbour site, and (ii) the possibility for tunnel- 
ing transitions to the nearest-neighbour site via a pair 
of virtual forward/backward moves. The first possibility 
contributes to the rate 7^, and the second one to the 
rate 7^. In case (i), the charge arrangements have two u 
charges on the path q and two v charges on the path q' , 
all of the same sign. The unrestricted integration over 
r in Eq. (|18|) yields all possibilities of mixing up the u 
with the v charges. The resulting sets T> x and T> x each 
have six different time-ordered arrangements as shown in 
Fig. H These sets make up the backward rate 7^2 an d 
forward rate 7^2? respectively. 



V x 

- - e e 

- e - e 

- e e - 
e - - e 
e - e - 
e e - - 



+ + © © 

+ © + © 

+ © © + 

© + + © 

© + © + 

© © + + 



A. Golden rule limit 

For very weak intersite coupling A s , the leading contri- 
bution in the series Q29[) has one u and one v charge. The 
respective ordered arrangements are + © and © + for 
the forward rate, and the corresponding charge-conjugate 
pairs for the backward rate. The noise integrals associ- 
ated have £1 = p\ = ±1 and read 



J x — i % S x 



dr e 



r 2K 



Equation J2jJ) with l(3"U|) for n = i = 1 then gives 
x 2 

7^ = e-^ [cos(wK) C x ± sinfjrif) S x } . 

With the substitution r = ±ix in Eq. Ij48(l . it is straight- 
forward to see that 



(48) 



(49) 



C x = tan(7riT) S x , 
and S x = cos(wK)T(l - 2K). Thus 7^ = 0, and 

j xx = eig sin(7rA") S x = 



2T(2K) 



(50) 



(51) 



This is the golden rule tunneling rate to the nearest- 
neighbour site in forward direction in the Schmid model. 
Clearly, the same expression 



it 



e Xg B sm(irK)S 1 



(52) 



holds for the rate in the SB model in the nonadiabatic 
limit. 



FIG. 1: The set T> x consisting of six different charge se- 
quences or paths constitutes the backward rate 7jT 2 - The 
charge-conjugate set "D 1 determines the forward rate 72V 



In case (ii), there are three charges on path q and one 
charge on path q' , and vice versa. Accordingly, there 
are three different sets denoted by £>2, 2?3, £>4, and the 
respective complex conjugate counterparts T>2, T>^, T>\. 
These six sets make up the rate 7^ (see Fig. 01. 



©2 

- - + © 

- - © + 

- © - + 
© - - + 



- + - © 

— h e - 

- © + - 
e - + - 



e 



+ - 

+ - © - 

+ e - - 

© + - - 



FIG. 2: The three sets T> 2 , T> 3 , X> 4 , each consisting of four dif- 
ferent arrangements, form, together with their complex con- 
jugate counterparts 2?3, the rate 7 1 ~ 2 - 

The quantum fluctuations are carried by three differ- 
ent arrangements of the cumulative charges p x , p 2 , p 3 , 
which are {1, 2, 1}, {1, 0, 1}, and {1, 0, —1}, and their 
charge conjugate counterparts. The corresponding con- 
nected complex noise integrals 



2,k 



^2,k i i& 



2,k > 



1, 2, 3 



(53) 
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in order A 4 , are found to read 



j± 



d 3 f e ±i ( T i +2r 2 +r 3 



'1'3 



T 2 T 12 r 23' r 13 



2 A' 



J+ = / d¥e ±i ^ +T ») 



'12 '23 



' 1 '2 '3 ' 13 



2 A 



2 A" 



Jf 3 - / d 3 fe ±! ( T i- T 3 



"/ x 2 A 






\ r l' r 3' r 12' r 23 / 


( T l r 3 ) 



2 A" 



The integrals J 2 2 and J+3 embody intermediate visits of 
diagonal states. The subtractions made serve to elimi- 
nate the reducible parts. The set T> 1 yields the backward 
rate "/ 22l which is zero. Upon transcribing the set T> 1 
into the {77, £} representation, we get the following linear 
relation between noise integrals J 2k (k = 1, 2, 3), 



ii-rrK t— 



^2,1 



J+ 
J 2.1 



-,-12-kK 



Jo 



+ e i2 * K J+ 2 



'2,3 



J+ = 



(54) 



The sets T> 2 , T> 3 , T> 4 yield complex partial "rates" 
7m P2I %A V ^ and 7i"ip4] - According to Eq. 
these partial rates vanish individually at T = 0. The 
resulting linear relations between the noise integrals are 



^ e -i2^A 



i4ttK\ t- 



J^,l 



i2ixK j- 1 t+ 
J 2,2 "T J 2,3 



j 2 : 2 + J+ 2 + (1 + e" 



^21^A^ 



0. 

J 2 - 3 = 0, (55) 



(1 



i2ttA n 



"^2,1 "t" "^2 



In addition, there are relations between the real and 
imaginary parts of the noise integrals J 2 x and J 2 2 . Upon 
rotating the integration paths by ±5, t = iix, the 
integrals J+j and J 22 receive a definite phase factor 
zfie^ i27rK . Thus we get 



C 21 = taa(2nK) S 2 ,i , 
C 2a = tan(2vr^) S 2j2 . 

Furthermore, for reasons of symmetry we have 



S. 



2,3 



. 



(56) 



(57) 



The set of linear equations l|54[) and (|55|l , together with 
expressions (|56*|) and (|57|) . yields the relations 



C 2 , 3 = tan(27rX) S 2<1 , 

COs(27rif) ] 5 2 ! . 



^2,2 = "[I 



(58) 



With the expressions l|56l) - l|58|l the noise integrals 
C2 1 j C2 2 1 C2 3 an d ^2 2 can be expressed in terms of the 
noise integral S 2 1 . As a result, also the individual partial 
rates in forward direction can be written in terms of the 
noise integral S 2 1 . The nearest-neighbour partial "rates" 



7l + 2 [P 2 +P 2 *] = - e ^sin(27r^)cos(27r^)5 2il , 
7it 2 [^3 + ^3*] = e^sin(27r^)[l + cos(27rX)]5 2il , 
7l + 2 [P 4 + ^] = -e^sm(2nK)S 2A . 

Now observe that these partial rates add up to zero, 

7 + 2 = . (59) 

Next, we recall that the set T> 1 in Fig. ^ determines 
the forward rate to the next-to-nearest-neighbour site in 
order A 4 ., 7^~ 2 , in the Schmid model. Using again the 
relations l|56l) - (|58|l . we obtain the compact expression 

7 2 + 2 - -e^ sin 2 (7rX) sm{2*K)S 2il . (60) 

Finally, we turn to the SB model. Because the signs 
of the u and v charges must alternate, the only charge 
arrangement contributing to the rate 7^ are the sets T> 3 
and f>l (cf. Fig. EJl. Thus we obtain 

x 4 

7+ =e-^sm(2nK)[l + Qos(2nK)]S 2A . (61) 

Upon combining the expression l|6UII with (|61|l we find 
the simple relation 



7 2 



1 SB 



\ 4 sin 2 (2nK) 7 + 2 . (62) 



v 2sin(7rlf)A s 
This concludes the discussion of order A 4 

C. The order A 6 



In the Schmid model, the irreducible clusters of order 
A 6 contribute to the rates 7+ , 7+ and 7+ . Now consider 
first the various noise integrals of this order. There are 10 
different complex conjugate pairs of possibilities to order 
the £ charges. Some of these pairs can be transferred into 
each other by reversal of the order of the charges. As a 
result, there are seven different arrangements left. The 
related noise integrals J 3k — C 3 k ± iS 3 k (k = 1, • • • , 7) 
are listed in Fig. The individual expressions can be 
be traced from Eq. (|31[) . The noise integrals JJ^ and 

J+2 are irreducible on principle, and we assume that the 
reducible parts of J+ 3 , • • • , J 37 are already subtracted. 
The cumulative charges of the noise integrals 



J 31 ,-- - , J 34 do not change sign. The real parts of 
these integrals are again directly related to the respec- 
tive imaginary parts. The corresponding relations are 
discovered upon rotating the integration paths by ±5, 
t = ±ix. The integrals then receive a definite phase 
factor ±ie Tl3 * K . We then find 



c 3,k = tan^Tr-FQ^fc i 



for k = 1, 2, 3, 4 . (63) 
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t± 
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FIG. 3: The different sets of cumulative charges, which (to- 
gether with the charge-conjugate counterparts) define the 
noise integrals Jg 1; • • • , Jg 7 . 



Now consider first the rate contribution 7 3 3 . The rel- 
evant charge sequences have three charges of type u and 
three charges of type v, all of the same sign. There are 20 
different possibilities to mix up the v charges with the u 
charges, i.e., there are 20 individual paths on the {q,q'} 
or {r],£} plane. Hence 7^3 emerges as a linear combina- 
tion of 20 noise integrals chosen from the list displayed in 
Fig. |3 and the coefficients can be traced from Eq. lj3TJ|) . 

Secondly, the rate 7^3 is made up of arrangements with 
four u charges (three of the same and one of opposite 
sign) and two v charges of same sign, and of the corre- 
sponding arrangements in which the u and v charges are 
interchanged. According to the four different possibilities 
to order the four u charges, there are four different com- 
plex partial "rates", denoted by 7^ 3 [fc] {k = 1, • • • ,4). 
There are 15 different possibilities to mix up the four u 
charges with the two v charges. Hence each of the par- 
tial rates 7^3 [k] can be written as a linear combination 
of 15 noise integrals taken from the set J^, • • • , J^ 7 . 
Taking into account also the complex conjugate arrange- 
ments, the rate 7^3 is made up by altogether 120 different 
charge sequences or tunneling paths. 

Thirdly, the rate 7^3 emerges from (i) the sets with 
five u charges (three of the one and two of the other 
sign), and one v charge, (ii) the corresponding complex 
conjugate sets, and (iii) the sets with three u and three 
v charges (each with two charges of the one sign and 
one charge of the other sign). In case (i) there are 10 
possibilities to order the u charges. This results in 10 
complex partial rates 7j t 3 [fc], k = 1, • • • , 10. Each is a 
linear combination of six noise integrals, according to the 
six possibilities to add in the v charge to the arrangement 
of u charges. Case (ii) yields the respective complex con- 
jugate counterparts. Case (iii) gives three possibilities to 
order the u and three possibilities to order the v charges. 
Accordingly, there are additional 9 partial rates 7^3 [k], 
k = 11, • • • , 19. Since there are 20 different possibilities 
to mix up the u with the v charges, each of these partial 
rates is a linear combination of 20 noise integrals, again 
taken from the list shown in Fig. In summary, there 
are 300 different charge sequences or tunneling paths al- 



together which determine the partial rate 7*3. 

All the noise integrals J 3 j (j = 1, • • • , 7) are related 
with each other such that 

7^3 = °. 

72 - 3 [fc] =0, k = l, •••,4, (64) 

7l " 3 [fc] =0, k = l, ■■■ ,19. 

With use of Eqs. 164(1 and (|63|l the noise integrals 
J^g , • • • , J^ 7 can be expressed in terms of the noise in- 
tegrals S 3 1 and S 3 2 - Upon using these relations the 
respective forward rates are found as 

7i + 3 = , 

72^3 - , 

x 6 (65) 
7+3 = e^sin 2 (7rif)sin 2 (27r!T) 

x [sin(nK)S 32 + sm(3irK)S 31 } . 

Finally, consider the forward rate in the SB model in 
order A 6 . Since the u and v charges must alternate 
in sign separately, there are only 32 different tunneling 
paths contributing to 7^, and the noise integrals involved 
are J^, J^g and J^ 7 - Expressing these in terms of the 
noise integrals S 3 1 and S 3 2 and using Eq. (|65|l . we find 

^ + -( 2sin t")A s ) 64sin2 ( 3 ^3- (66) 

In conclusion, we have found again formidable cancel- 
lations between the contributions of the various tunneling 
paths in the Schmid model. All those rates of order A| 
which involve pairs of virtual forward/backward moves 
vanish. In addition, we found again a simple relation be- 
tween the rate in the SB model and the one nonvanishing 
rate in the Schmid model of same order in A. 



D. Higher orders 

With increasing order of the perturbative expan- 
sion, the number of contributing tunneling paths in- 
creases enormously. Therefore, the computation, al- 
though straightforward, becomes rather elongated. The 
general strategy is based on two observations: 

(i) The set of equations l|47|) results in a set of linear 
relations between the noise integrals jf~ (£ ) of same order 
I and different £. 

(ii) Consider an arrangement of "fluctuation" charges 
£ = ^ ( d ) f or w hich the bias phase <p t has a definite sign. 
We see from Eq. (|32f) that the sign of ip e is definite when 
all cumulative charges Pj have the same sign (apart from 
zeros), say sgn^-) = +. It is now decisive that the 
real and imaginary part of each individual noise integral 
J^{£, ^ ) are related to each other. To determine the re- 
lation, we rotate the integration paths by zt^, f — ±ix. 
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With these rotations, ^ ) picks up a definite phase 

factor, which is q=z e ±llr ' 1_Ar ^. Thus we find 



•(d) 



t(m(£irK)S e (ti 



"(d) 



(67) 



This is the generalization of Eq. (|63l) to order A . 

Property (i) is in general a virtue of the Ohmic scaling 
limit, Eq. It also holds at finite temperatures, where 
Eq. I|47(l is replaced by Eq. I|20[l (see remark at the end 
of Section ITTT|l . Property (ii) is based on the scale invari- 
ant logarithmic charge interaction H46JI . It does not hold 
when scale invariance is broken by finite temperature. 

Upon combining (i) with (ii) we find for fixed I that 
all those noise integrals J'f 1 (£ ) , which have one or sev- 
eral of the cumulative charges pj equal to zero, can 

be expressed in terms of noise integrals S e (^ ) in which 
all Pj are nonzero, which entails that the respective p. 
have all the same sign. Making use of these relations, 
we discover formidable cancellations among the various 
path contributions to in Eq. ij^UI) . In fact, all 

tunneling paths which describe detours via virtual for- 
ward/backward hops cancel each other. Thus we find 



<t = o 



for £ > 



(68) 



Only direct tunneling paths which have all u and v 
charges positive contribute. Finally, we obtain in gen- 
eralization of Eq. (|65|l 



In — ln,n — e ^"^r^n ' 



U+ = 



(-iy 



2n-l 



2n 



(69) 



{«*}' 



fe=i 



Next we turn to the discussion of the SB model. Mak- 
ing again use of the relations between the noise integrals 
of same order n, we discover that the weight of the 
SB partial rate 7+ = ex^W^ given in Eq. (|4"5|) can di- 
rectly be expressed in terms of the weig ht U+ n of the 



rate 7+ = ex^ n U^ n in the Schmid model. We find 



( ' [2sin(7riq] 2 " n - r 



(70) 



For later convenience, we now introduce a frequency 
scale e analogous to the Kondo scale in the Kondo 
model and the scale T B in QlPsiiS The relation of e 
to the bare parameters of the Schmid model is estab- 
lished by the exact self-duality of this modeli As far as 
the SB is concerned, it is convenient to absorb the factor 
[2 sin(irK)]~ 2n occurring in Eq. (|70|) into a corresponding 
redefinition of the Kondo scale. Thus we have 



2-2if 



2 2 - 2K tt 2 A 2 



T 2 (l-K) A 2 



SB 



T 2 (K) 



,2K 



2 2K 



,2 if 



(71) 



At fixed renormalized coupling e , the SB bare coupling 
A SB is then related to the bare coupling A s of the Schmid 



model as 



A 



SB 



(72) 



Upon employing relation l|7()(l . we then find that the 
partial rate 7+ of the SB model is simply given by 

7+ = (-l)"- 1 4sin 2 (n 7 r^)7+. (73) 

The intriguing functional expression (|T3I) is a major result 
of this paper. It relates the incoherent partial tunneling 
rate of order A|g in the SB model to the full tunneling 
rate from site to site n in the Schmid model. 

At the end of this subsection, let us briefly pause to 
shed light on the relations (|71|l and (|72|l from different 
perspective. First, we note that the relation l)72|l is in 
agreement with findings from a unification of the SB and 
Schmid model within the quantum group SU(2) q m& Sec- 
ondly, we recall the relationship of the Ohmic SB model 
with the Kondo model. While the original isotropic 
Kondo model is at K — 1, the anisotropic Kondo model 
with independent exchange constants pj^ and pj ± de- 
scribes variable K. The correspondence between the 
two models is universal, i.e. it is independent of the 
regularization prescription chosen for the singular scat- 
tering potential^ in the regime pj ± = A gB /w c 1 
and p\J\\\ = |1 — K\ <C 1. In the parameter range 
pj ± pj|| <C 1, the Kondo energy (modulo a factor 
2 /it) is expressed in terms of the bare parameters asi^ 



pJ± 
2p,h 



1 



A 



SB 



(1 - K) (2w c ) 



A" 



(74) 



The second form ensues from use of the correspondence 
relations. Now observe that the second form in Eq. I|71() 
coincides with the expression (|74|l in the regime around 
the critical coupling, |1 — K\ <C 1. These two independent 
conclusions strongly confirm the above analysis of the 
perturbative expansions of the Schmid and SB model. 



E. Statistical fluctuations 

From the characteristic function (0 at T = we see 
that the connected moments of the probability current 
( / > c = E„ $ = £„ n lt take the form 



r( m )\ 



(75) 



where what is called (/' m ^) c corresponds to (n m (t)) c /t. 
Now, since l£/e is of order Ag n e( 2K-2 )", as follows from 
Eq. H69|l with Eq. I|28(l. simple moment relations can be 
derived. We find 



(J(r. 



(I)c 



a s d 

2 OA, 



2(K-l)de) e ^ c 



gyn-l 1 



(76) 
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Thus all statistical fluctuations of the transport process 
at T = can be extracted directly from the current. 

We conclude with the remark that these findings from 
the nonequilibrium Keldsysh approach are in agreement 
with results from the thermodynamic Bethe ansatz ^A 



VI. FROM WEAK TO STRONG TUNNELING 

A. BSG or Schmid model 

There is an exact duality for the nonlinear mobility at 
temperature T of a Brownian particle in a tilted periodic 
potentialiA^ In the duality transformation, weak and 
strong corrugation of the potential are exchanged and the 
spectral density is transformed. One finds that super- 
Ohmic damping maps on sub-Ohmic and vice versaji*i 
while Ohmic damping maps on Ohmic (with different 
high- frequency cut-off), but the damping strength K 
maps on 1/ifiii The duality becomes an exact self- 
duality in the Ohmic scaling limit A Independently, ther- 
modynamic Bethe ansatz computationsAS established 
that there is an exact nonperturbative self-duality for 
the current in the BSG model at zero temperature. Con- 
sidering the dimensionless current T(K, e/e ) = 2tt Kl/e 
as an analytic function of K, gives the self-duality as 



X(K,e/e ) = l-l(l/K,e/e ) 



(77) 



This equation applies for all K > 0. Doing perturbation 
theory in the strong-backscattering limit yields 



1(K, e/e Q ) = Y^i n {K) 



(2K-2)n 



(78) 



The thermodynamic Bethe ansatz gives the coefficients 
i n (K) in the formAS 



,(*) 



(_i)n-i r(|)r(i + nK) 



T[l + n(K-l)[ 



(79) 



With the perturbative series (J7SJ), the self-duality l(77|) 
then yields the corresponding asymptotic or weak- 
backscattering expansion. 

The existence of the self-duality Ij77(l has prompted 
questions about the analytic structure of the BSG trans- 
port problem. It has been argued that, knowing a priori 
that Eq. I|77(l holds, it should be possible to find the exact 
expression for X without using the Bethe ansatz. Follow- 
ing these lines, exact integral representations for the cur- 
rent were derived which harbour both the perturbative 
and the asymptotic expansion. 1 - 15 - 16 

Using the explicit form of the rapidity-dependent 
transmission coefficient in the thermodynamic Bethe 
ansatz, the shot noised and also the full probability distri- 
bution P(fc, t) were derived^A The strong-backscattering 
expression for the characteristic function was found to be 



in Poissonian form 



P(k,t) = Y[cxp[t(e lkr ' 



1)7+ 



7.,; 



r(§)r(*fn) 



n! r[| + [K - l)n] 2tt \e 



(80) 

(2K-1)n 



While a classical Poisson process requires all the 7+ to be 
positive, the signs of the rates are alternating. The first 
contribution is indeed a Poisson process for the tunneling 
of electrons, but the tunneling process of a correlated pair 
of electrons (and multiples thereof) has the wrong sign, 
which is a distinctive feature of quantum interference. 
The weak-backscattering expression was found as 



P(k,t) 



In 



Jktt/2-irK 



J] expke- ifc "/*_ 1)7+1 , (81) 



n=l 



(-1)"- 1 r(|)r(f) 

nl r[| + (i-l)n]27r/v Ve 



(*-3)n 



The first exponential factor represents the current in the 
absence of tunneling. The exponents have now the term 
e -ikn/K _ rp^g £ ac ^ or \JK = v is the signature of tun- 
neling of Laughlin quasiparticles and multiples thereof. 
The minus sign occurs because their tunneling diminishes 
the current, instead of building it up as in the strong- 
backscattering limit. The sign of the transfer weights 7+ 
is now dependent on K. Using the reflection formula for 
gamma functions, the sign of 7+ is the one of cos(nir / K) 
and therefore positive for modest n when K ^> 1. As 
K — > 00, the classical regime is reached, in which all 
rates 7+ are positive, and the fluctuations disappear. 11 
In this limit, the expression (|81|l is summed to the form 



P(*,t)=«p|ifct 5 ^^l-(^) |6(6-e ). 

We see from Eq. (I80|) that 7+ is exactly of order Xg n . 
Hence the integrable approach is in correspondence with 
the findings fflfy of the perturbative analysis. Conversely, 
the rigorous nonequilibrium Kcldysh approach presented 
in Sec. thus corroborates the thermodynamic Bethe 
ansatz. As an important by-product, Eq. (|80|l provides 
an analytic expression for the particular linear combina- 
tion of noise integrals given in Eq. I|tj9|) • 

We have mentioned already at the beginning of Sec. El 
that the noise integrals l7^(£) diverge at short distances 
for K > 1. At this point two remarks are appropriate, 
(i) It is most natural to define the regularized integrals 
for if > i as the analytic continuation of their values 
for K < i. Since the integrals have simple poles and 
there are no branch points, the analytic continuation is 
well-defined, (ii) The friction factors exhibit zeros at the 
locations of the poles of the integrals. As a result, the ex- 
pression I|tj9fl , in which friction factors and noise integrals 
are combined, is a smooth function of K for all K, as can 
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be seen from expression (|80ll . Also the Kondo scale for 
the BSG model, Eq. (|71|) . is nonsingular for all K. On 
the other hand, the remaining singularities in the SB ex- 
pression H7U|) can be fully absorbed into the definition of 
the Kondo scale for the SB model given in Eq. JTlJ, as 
we see from relation (J2HJ). From this we may conclude 
that the rate expressions found from our analysis for the 
Schmid and SB model are valid for all K > 0. 



B. Spin-boson model 

With use of relation l|73|l and expression (|80|l for the 
BSG rate 7+, we immediately obtain an exact expression 
in closed analytic form for the partial tunneling rate 7+ 
in the spin-boson model. With this, the weak-tunneling 
expansion of the full rate 7 + , Eq. (|44|) . is found to read 



(2K-2)m 



1 T{Km) [1 - cos(2nKm)] 
ml T[l + {K-l)m] 



(82) 



For rational K, the series (|82|) can be written as a linear 
combination of hypergeometric functions. 

In the regime K < 1, the perturbative series l|82|l is 
absolutely converging for large enough bias. For K > 1, 
the series is absolutely converging for small enough bias. 
The leading term is the golden rule rateia 



the order of integration and summation. This yields 



= Ue-1,/'^ 

2m 



z — 1 — z K u 2 — \jz — 1 — z K u 1 



(86) 



The integral representation (|86|l converges for all K 
and all e/e , not just in the regime where the series l|82|) 
does. Evidently, Eq. I|86(l is generally valid. Inter alia, 
it must also hold the corresponding strong-tunneling or 
asymptotic expansion for small (large) enough bias when 
K < 1 [K > 1). 

Consider first the asymptotic expansion for K < 1. 
Upon changing variable z to y — z 1 ~ K fu x and y = 



,1-K 



/u 2 , respectively, Eq. (|8u[l transforms into 



7+ 



1 



~o_ [ 

1-K 2m Jc y 



x Re 



' y - 1 - [ye 



(87) 



l-y~> 



Next we expand the curly bracket for small e/e and in- 
terchange summation and integration. Using Eq. I|84|l . 
the strong-tunneling series is found as 



7gr 



2T(2K) 



SB 



sm 2 (irK)T(K) ^ 1 e 



V^ra + k) 



2K-1 



(83) 



?;.=0 



0+ — 



2n 



For K < |, the coefficient b n {K) is given by 



The coefficients a m (K = -|) are zero when m/1. Hence 
7gr coincides with the full rate for K — \. 

An integral representation for 7+ is found by writing 
the quotient of the Gamma functions in Eq. I|82l) as a 
contour integral^ 

(84) 

T{x + y) 2m J c 

The contour starts at the origin, circles around the 
branch point in counter-clockwise sense, and returns to 
the origin. Equating x with Km and y with | — m, and 
putting u 1 = (e/eo) 2K ~ 2 and u 2 = e z2 ' KK u 1 , we get 



7 T = 



v m= 1 



r(m-|) 



Re 



ml 



2m Jc z V z — 1 



(85) 



Since the weak-tunneling series l|82|) is absolutely conver- 
gent within the circle of convergence, we can interchange 



b n {K) = 2^[^ R {\-n)]c n {K) 



where 



c n (K) 



m-n 



- JL -} 



n\ (I-njrKl-n)^] 



(89) 



(90) 



For I < K < 1, the contour does not circle the branch 
cut of the first term in the curly bracket in Eq. I|87(l . 
Hence only the second term contributes and yields 



b n (K)=c n (K) 



(91) 



The strong-tunneling or asymptotic series (j88(l is an 
expansion in powers of e 2 /eo- The leading term repre- 
sents the forward rate of the symmetric SB model, 



7o + 



A K ristrby] 



(92) 



2(1-K)\ 



where A K = 2 sin 2 [kK/2(1-K)} for K < ± and A K = 1 



for i < K < 1. 
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K = 2/3 




FIG. 4: The scaled rate 7 + /7qr is plotted versus Xg B for 
various K < 1. The circles represent the weak-tunneling se- 
ries and the squares the strong-tunneling expansion. The full 
curve is the hypergeometric function expression 




FIG. 5: The same scaled rate, but now for various K > 1. 



Consider next the weak-damping limit. For K « 1, 
the perturbative series ((821) as well as the asymptotic 
series l|88f) can be summed to the simple form 



is useful to divide the damping regime into the sections 
K = p + k with p = 1, 2, • • • and < k < 1. In the first 
step, we change in the integral representation ((86(1 from 



variable z to variable t — 



U 1/2 Z 



K 



. This gives 



7 + 



Re 



1 e 
~K 2tTi 



At 



H{t, Ul ) 



(94) 



H{t, Ul ) = Jt- l + e-a+^/^v) 1 /^ 



(95) 



Second, we expand H(t 1 u 1 for small u 1 and interchange 
integration and summation. Using 1)84(1 . we then obtain 
the asymptotic series as 



7+ 



2^F E dm(K) 

v m— 1 



(2/K-2)m 



(96) 



d m (K) 



(-l) m 2T(f )sin[i±2m7r]sin(f tott) 
~\ KY\\ + (j{ - l)m] 



The asymptotic series (jMJl for the SB model in the regime 
K > 1 bears a strong resemblance with the corresponding 
series of the self-dual Schmid model (|81|l . The powers of 
e/e follow from those of the perturbative expansion 1(82(1 
by the substitution K — > 1/if, as also do major parts 
of the coefficient d m (K). However, complete duality is 
spoilt by the alternating sign and by the sine-factors. 

Fig. [3 shows plots of the normalized rate 7 + /7(t R for 
K > 1. The normalized rate goes through a maximum 
which is shifted to higher x SB when K is increased. At 
large enough x SB , the rate 7" 1 " falls below the golden rule 
rate. Hence there is constructive interference of tunneling 



at small and intermediate 



J SB' 



and destructive interfer- 



ence in the large tunneling regime. 



j+=nKA 2 SBi 



(93) 



which is in agreement with findings from a systematic 
weak-damping self-energy methodic 

In Fig. 01 the normalized rate 7 + /7q R is shown for 
different K < 1. The horizontal line represents the par- 
ticular case K = |. For K < i the full rate is always 
lower than the golden rule rate. Hence the numerous tun- 
neling paths interfere destructively in this regime. For 
i < K < 1, the tunneling contributions interfere con- 
structively for all a; SB so that the full rate is always above 
the golden rule rate. In the regime 4 < K < |, the nor- 
malized rate goes through a maximum as tunneling is 
increased, and finally falls below one. This reflects con- 
structive interference at small and intermediate x SB , and 
destructive interference at large x SB . 

Next consider large damping K > 1 and large bias. It 



C. Conjecture on decoherence at T — 

So far, we have left aside the damping of the coher- 
ent oscillations of the populations. If one aims at calcu- 
lating frequency and decoherence rate of the oscillation 
directly, the A-dependence of the irreducible clusters dis- 
cussed above has to be taken into account. The calcu- 
lation of the respective self-energy and of the decoher- 
ence rate has been carried out in special limitsji e.g. for 
K <C l. 18 Here we pose the conjecture that the strong 
tunneling expansion of the tunneling relaxation rate 7 + , 
Eq. ((88(1 . and the corresponding expansion of the deco- 
herence rate 



v( d ) 



^ ' 

n=0 



v( d ) 



2V^ 



fn(K) 



(97) 
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are closely related, viz. that the ratio of the coefficients 
of these expansions is given by 



fjK) 



sm 



2{l-K) 1 



<K< 



Equation l|98|) is in agreement with all exact expressions 
known in special limits. 

The decoherence rate of the unbiased SB model 
has been calculated within the framework of integrable 
QPT 19 The respective result (translated into our nota- 
tion) is 



(d) 

7o 



i sin 



tK 



2(1-K) 



£[ 



K 



2(1- K)i 
2(1-K)\ 



(99) 



Corresponding expressions hold for the current (I) c and 
for the moments (/(™ l ') c given in Eq. 1)76(1 both in the 
weak- and strong tunneling regime. Equally, the thermal 
enhancement of the tunneling rate and the decoherence 
rate in the SB model has a form analogous to Eq. (|101fl . 

We see that the leading enhancement at T < £ varies 
with T 2 . The power 2 is a signature of Ohmic dissipation. 
In addition the prefactor is universally given by the curly 
bracket with the second order differential expression in 
Eq. (|1()1J) . The universal form is closely related to the 
Wilson ratio in the AKM, in which the specific heat is 
related to the static susceptibility^ The physical origin 
of the T 2 enhancement is the low- frequency thermal noise 
of the dissipation^ 



is in correspondence with expressions 1)98(1 and (|92[1 for 
K < |. The relation 198fl is known to hold also in the 
biased case in the regimes K <C 1 and K close to 

These checks suggest that the conjecture (|98|l holds 
for the SB model at zero temperature in the regime < 
K < |. Upon combining Eqs. ® - HUB with (JHSJ the 
weak-bias expansion of the decoherence rate is found as 



20? 



n=0 



[ irK I 1 
ll-K\2 



n)]c n {K) 



(100) 



If the conjecture (|98|l is correct, then the series expression 
(|100fl represents the lower bound for decoherence in the 
SB model in the scaling limit in the regime < K < i. 



D. Leading enhancement at low temperatures 

We now extend the discussion to the asymptotic low- 
temperature regime. The leading thermal correction in 
the noise amplitude factor G t is found from Eq. 127fl as 

2 2l ~ X 2 

G e = Gtf = 0){l-Jfy^ 2 (^ Pi T 3 ) + 0(tf 4 )} , 

i=i 

The important point now is to recognize that the an- 
noying term in the round bracket is the bias phase (|32[l . 
Therefore, this term can be generated by differentiation 
of the noise integral (|3*T1) with respect to the biasi^S We 
thus have, for example, 



lt{^T)= 1 



7T 2 „ 

K-T 2 



d 2 



(def 



T +(e,T = 0). (101) 



VII. CONCLUSION 



In summary, we uncovered exact functional relations 
between perturbative integrals occuring in the nonequi- 
librium Keldysh approach applied to the Schmid and the 
SB model at T = 0. With that and with use of the corre- 
spondence of the Schmid with the BSG model, we found 
agreement between the results of the Keldysh approach 
and the thermodynamic Bethe ansatz for the BSG model. 
As a second result, we discovered exact functional rela- 
tions between the rates in the Schmid model and the 
partial rates in the SB model. Starting out from the 
perturbative series obtained by use of the functional re- 
lations, we put up an exact integral representation for 
the tunneling rate in the SB model, from where we de- 
rived the asymptotic or strong-backscattering expansion. 
Thirdly, we conjectured and checked a relation between 
the relaxation and the decoherence rate in the SB model. 
Finally, we calculated the leading thermal enhancement 
of the transport and of all statistical fluctuations and 
found universal behaviour. 

Our results are relevant, e.g., to charge transfer, 
macroscopic quantum coherence (qubits), and diverse 
quantum impurity problems. 
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